High-order noise filtering in nontrivial quantum logic gates 
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Treating the effects of a time-dependent classical dephasing environment during quantum logic 
operations poses a theoretical challenge, as the application of non-commuting control operations 
gives rise to both dephasing and depolarization errors that must be accounted for in order to 
understand total average error rates. We develop a treatment based on effective Hamiltonian theory 
that allows us to efhciently model the effect of classical noise on nontrivial single-bit quantum logic 
operations composed of arbitrary control sequences. We present a general method to calculate the 
ensemble-averaged entanglement fidelity to arbitrary order in terms of noise filter functions, and 
provide explicit expressions to fourth order in the noise strength. In the weak noise limit we derive 
explicit filter functions for a broad class of piecewise-constant control sequences, and use them to 
study the performance of dynamically corrected gates, yielding good agreement with brute-force 
numerics. 



Dynamical error suppression strategies have been 
demonstrated as a means by which errors due to deco- 
herence may be suppressed during qubit memory opera- 
tions [THS] . In such filter-design framework [51115] 
has successfully shown how to precisely estimate average 
error rates in the presence of classical, time-dependent 
noise. Expanding this analysis beyond the identity oper- 
ator has proved challenging due to the need for efficient 
techniques to treat a random, time-varying noise term 
that does not commute with the applied control field. 
Understanding the influence of such time-dependent pro- 
cesses during control operations is vital, however, as en- 
vironmental decoherence sets the lower-bound on achiev- 
able gate error rates in a quantum informatic setting. 

In this letter we address the challenge of character- 
izing and mitigating decoherence due to classical noise 
during nontrivial single-qubit operations. We calculate 
the ensemble-averaged entanglement fidelity for an ar- 
bitrary control sequence to fourth order in the noise 
strength, incorporating terms to the third order of the 
Magnus Expansion in an effective Hamiltonian treat- 
ment. For concreteness we explicitly calculate the filter 
functions to lowest nontrivial order for sequences com- 
posed of TT-rotations about Cartesian axes with arbitrary 
rotation rates - a class including dynamically corrected 
gates (DCGs). Our results permit detailed calculation of 
the generic (amplitude and phase) errors that result from 
applying a quantum control operation in the presence of 
pure-dephasing noise, and validate perturbative predic- 
tions of an increased order of error suppression [141 115j . 
This straightforward analytical approach is compared 
against brute-force numerical calculations of the evolu- 
tion of the Bloch vector, and shows excellent agreement. 

We consider the canonical dephasing environment in 
which the system evolves freely under a Hamiltonian of 
the form H = }^\^l-\-rj{t)]az, where $7 is the unperturbed 
qubit splitting, 77 is a time-dependent classical random 
variable, and is a Pauli operator. In the case of free- 



evolution the presence of a nonzero ri{t) produces de- 
phasing in an ensemble average. However, during driven 
operations where one applies a control field proportional 
to Gx or Gy, the presence of a pure-dephasing noise envi- 
ronment yields both polarization damping and dephasing 
effects. Both must be considered in a full treatment of 
gate-errors. This accounts for the most significant cor- 
rectable forms of decoherence in experiment; most re- 
maining polarization-damping errors are due to stochas- 
tic processes that cannot be corrected through dynamical 
error suppression. 

We begin with an outline of our method. The to- 
tal Hamiltonian (in the rotating frame at Q) is H{t) — 
HQ{t) + Hcit), over t G [0,r]. The operator Ho{t) = 
'7(^)cz/2 represents a time-varying dephasing environ- 
ment, while Hc{t) describes an interaction between the 
system and an external control device that, in princi- 
ple, may be used to implement arbitrary rotations of the 
Bloch vector. In general, evaluating the total propagator 
U{t) = Texp{-i f^H{t')dt') explicitly for iJo(0 7^ is 
difficult due to noncommuting terms in H{t). We there- 
fore proceed by factoring out that part of the qubit evo- 
lution that is due solely to the control and expressing 
the residual 'error propagator' U (r) in terms of a time- 
independent effective Hamiltonian that can then be eval- 
uated, following a general procedure laid out in Ref. [2]. 

Defining the propagator Uc{t) ~ 
T exp{—i Hc{t')dt'), it can be shown that 
U{t) — Ul{t)U{t) satisfies the equation of motion 
.d|^t) ^ Ho{t)U{t), where H^{t) = Ul{t)Ho{t)U,{t) 
[ini [H]. If Hc{t) enacts a target unitary operation Q 
in the absence of decoherence, then we may write the 
total noise-affected operation as U{t) — QU{t). The 
error propagator can then be expressed in terms of a 
time-independent effective Hamiltonian H, defined by 
C/(t) = cxp{—iHT). 

For a dephasing Hamiltonian we have Ho(t) = 
'^Ul{t)azUc{t)\ the term Ul{t)azU^{t) is a rotation of 
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the Oz operator due to the control. We may therefore 
define a time-dependent 'control vector' in a Cartesian 
basis {I = x,y,z), si{t) = ^21^1,1^^ where |si(t)| = i, 
yt e [0,r], such that Hoit) = ^si{t) ■ a. Since Hoit) 
belongs to the Lie algebra sm(2), Vi € [0,t], the effective 
Hamiltonian H derived from it is also in su{2). 

We may now write U{t) = (5exp[— ia(T) • (?], where 
the effect of the noise on the ideal operation Q is en- 
coded in the 'error vector' a(T). Here a(T) — ai{T)l 
represents error contributions to Q due to terms in H 
proportional to the Cartesian components of the Pauli 
operator. In general, the action of the control converts 
pure-dephasing noise to arbitrary rotations of the qubit's 
Bloch vector, producing both dephasing and polarization 
damping errors. 

In most cases, an explicit expression for exp[— m(T) • a] 
can only be found using approximation methods. For 
our purposes, it is desirable that the simple exponential 
form of the operator be retained to any level of approx- 
imation. This can be achieved through the use of the 
Magnus expansion [TSl[Tn]. Using this method, the ex- 
ponent is expanded in an infinite series of time-integrals 
over nested commutators of Ho{t) at different times. In 
conjunction with the identity [u ■ a,v ■ a] = 2i{u x v) ■ a, 
u, i/ e M'^, we find that we can write a(T) • a = J^'iLi '^i ' ^ 
to all orders i. The first three terms in series expan- 
sion of the error vector are ai(T) — ^ dtr]{t)si{t), 

^2{t) = zJQdt2 J^^ dtiri{ti)ri{t2)s2{ti,t2), and a3(r) = 
hlo '^^3 Jq^ dt2 Jq'' dtir]{ti)r]{t2)'n{t3)s3{ti,t2,t3). Here, 
high-order commutators in the Magnus expansion have 
been reduced to vector cross products of the con- 
trol vector at different times, 5*2(^1,^2) = 5*1(^2) x 
siih) and 53(^1, ^2,^3) = Slits) x [31(^2) x si(ii)] + 

[51(^3) X 8*1(^2)] X Si(ti). 

We evaluate the net effect of the gate operation Q 
in the presence of noise via the ensemble average en- 
tanglement fidelity [10], {T{t)) = <^|Tr (QJ7(r)/2)|^^ 
which, when written in terms of the error vector, be- 
comes (J"(t)) ^ \ [(cos[2|a|]) -1-1]. To evaluate the fi- 
delity, we write \a\ — a^)^/^, express the cosine term 
as a Taylor series and substitute the Magnus expansion 
ai = ^Y^=\ '^or each component of the error vector. 
The result is an infinite series of multi-dimensional in- 
tegrals over products of multiple-time noise correlation 
functions and components of the control vector si{t). For 
example, the lowest order effects of the noise are captured 



by (all) = r^IZ'iAti)siAt2){vitiHt2))dhdt2, 

for I = x,y, z. 

Assuming the noise is Gaussian, only correlation func- 
tions {r]{ti)...ri{tn)) for which n is even contribute. Fur- 
ther, applying the Gaussian moment theorem, each of 
these can be written in terms of simple two-point corre- 
lation functions. Using the root mean square deviation 
of the noise, A77 = (rjit)^), as a measure of the noise 
strength, we can define a parameter ^ = A77T/2 which 
provides an upper bound for the magnitude of each term 
in the series expansion of the cosine function [51] [55] . If 
we restrict our analysis to weak-noise/ efficient-control 
conditions under which ^ < 1, then higher-order terms 
provide diminishing contributions to the total error and 
we may truncate the series. To fourth order in ^ we find 
that 

{cos[2\a\]) = l-2e{{dh)} 

(3(4,) + 2(51,^3,)) - {ahali,)^ (1) 

where sums are implicitly performed over I and I' , I I' , 
and the tilde indicates that the magnitude of the integrals 
has been factored out (—1 < dij < 1). This expression 
includes terms to third order in the Magnus expansion 
as they contribute to the same order in ^ as the second 
order term. 

The various contributions to Eq. [T] may be calculated 
explicitly by Fourier transforming the noise and control. 
For instance, we may write 

(aL> = ^ ri?/i,/Hr^d- (2) 
47r 7o 

with yij{uj) = —iijj si^i{t)e^'^*^dt, capturing terms 
proportional to uj. We sum over / to write i^i(w) = 
— 12i\yi,ii^)\ J corresponding to the terms in 
the first line of Eq. [T] Similarly, we account for the pro- 
liferation of higher order terms arising from the vector 
cross product by defining i^p 2(w, w', r), where p is an 
index over terms proportional to in Eq. [l] These 
terms contain four-point correlation functions in time, 
{ri{ti)ri{t2)ri{ts)ri{t4)) ^ which may be explicitly evaluated 
using the Gaussian moment theorem. 

We may then compactly write the entanglement fi- 
delity 
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Terms to all orders in ^ may be evaluated using the same procedure. 
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With these expressions we have reduced the effect of 
a time- dependent dephasing environment and a time- 
dependent control Hamiltonian to integrals incorporating 
only stationary statistical properties of the system: the 
noise power spectral density and the spectral functions 
to arbitrary order, Fi. These terms contain all relevant 
information about the applied control, and the explicit 
inclusion of terms proportional to all ai captures both de- 
phasing and polarization-damping errors produced dur- 
ing control operations. We refer to the terms as the 
filter functions for the total control operation, in anal- 
ogy with previous work on dynamical decoupling [TUHT^ . 
The leading nontrivial term is closely related to the idea 
of spectral overlap functions between control and noise 
that has been studied previously [23Hli] , but a gener- 
alized derivation to high order has not appeared to the 
best of our knowledge. 

We now consider a specific case for concreteness that is 
germane to many coherent control experiments, includ- 
ing the implementation of dynamically corrected gates. 
We define Hc{t) as piecewise-constant such that, during 
the j-th interval, the control is intended to execute ai^ , 
restricted here to either the identity /, or a rotation of 
the qubit Bloch vector through ±7r about one of the three 
Cartesian directions. In this notation, Ij = /, x, y, z. 

The control propagator may be written explicitly such 
that during the j'th driven operation we have Uc — 



exp[ 



tj^ijai- /2\a\fj^i, where fi^'' gives the 
driven rotation rate about axis I in time-bin j. The op- 
erator cTvj-i = <Jij_i(yi-_2---Cfi^ describes the cumulative 
effect of all completed rotations in the preceding time 
segments. 

We restrict our presentation to terms in the entan- 
glement fidelity of order appropriate for the case of 
weak dephasing noise. Higher order contributions are 
straightforward to calculate explicitly, but involve many 
dozens of terms. In this case we only require first-order 
components of the error vector and may approximate 

(J"(t)) « \ e"2^<"i.'> + 1 . Using the first-order filter 
functions we have {J-{t)) « ^ [exp(— x(r)) + 1], where 

For piecewise-constant control as described above, the 
terms contributing to Fi^i{uj) may be written: 



k 



.2 - / n^i^ 



j2 - f7 



(4) 



where Sap = 1 ior a — f3 and is zero otherwise. Here N^'^ 
represents the number of times ctq or appear in the 
sequence up to the j-th interval. We note that the prefac- 
tors in these equations are similar to those derived from 
the steady-state master-equation treatment of a driven 
two-level system in the presence of dissipation |26| . 

Using this approach we are able to analyze the ef- 
fects of classical noise on a complete set of both prim- 
itive (standard) and dynamically protected single-qubit 
gates of interest for quantum logic. As an example, in 
Fig. [l] we show the total filter functions derived from this 
treatment for X operations in primitive and DCG formu- 
lations [m [57]. Consistent with previous perturbative 
treatments of quantum-mechanical baths, we find that 
the order of error suppression, given by the low-frequency 
roUoff of is increased in the DCG relative to the 

primitive gate [131 HH] . The extension of the duration of 
Xdcg relative to X is manifested as a decrease in wpi, 
the frequency above which the filter function takes value 



unity and noise is passed largely unimpeded. By exam- 
ining the phase and amplitude components of the filter 
function independently (i^i,z and Fi^x), we see that for a 
DCG we improve the order of error suppression primar- 
ily in the term that commutes with the control operation 
(e.g. amplitude for an X gate). 

The form of the filter functions for other operations 
and their DCG constructions (e.g. Zg, H) show similar 
behavior. Dynamical decoupling (dynamically protected 
/) may also be treated using Eq. Q, and will be ad- 
dressed in detail in a separate manuscript. 

Using the lowest-order approximation for the entangle- 
ment fidelity and the specific forms of the filter functions 
for Xbcg we calculate the error probability for differ- 
ent noise environments and plot these in Fig. [2j As ex- 
pected, in an environment given by S{uj) = a/uP', with 
a a scaling factor (c.f. Ref (Mj): we find significant ben- 
efits from using the DCG construction. By contrast, in 
a white noise environment with a sharp high-frequency 
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cutoff {S{u!) = aQ{u!c—u!)) the significant high-frequency 
spectral components of the noise and the extended dura- 
tion of the DCG construction can yield a net performance 
degradation in the event of long and large ujc- The 
relationship between these two quantities thus provides 
a dominant practical limit on the applicability of DCG 
construction in realistic settings; so long as < the 
DCG provides net performance enhancement. 

Validation for the assumptions and approximations 
made in this treatment comes from performing detailed 
brute-force numerical simulation of the evolution of the 
Bloch vector in the presence of a noisy environment, and 
averaging over multiple trajectories, r]{t) (Fig. [2];). The 
trajectories are chosen to exhibit the statistical proper- 
ties of a desired power spectral density. Calculations in- 
volving only the terms in Eq. [4] accurately reproduce nu- 
merically calculated error probabilities to within ^ 20% 
for the X gate. In complex DCG sequences efficient de- 
coupling reduces leading-order terms until residual con- 
tributions become comparable to higher-order terms. We 
observe that the lowest-order filter functions underesti- 
mate error and deviation from numerics grows with more 
complex sequences or stronger noise, but remains a factor 
of order unity for the cases we have studied. 

In summary we have developed a theoretical treatment 
permitting the calculation of error rates due to time- 
varying classical noise during arbitrary control opera- 
tions, to arbitrary order. We have explicitly produced 
a general high-order approximation to the qubit's entan- 
glement fidelity and using the first three orders of the 
Magnus expansion. Building on these results we have 
presented an example case of piecewise-constant control, 
and given simple, leading-order expressions for the fil- 
ter function. This allows any experimentalist to evaluate 
expected error rates to the correct order of magnitude 
during nontrivial control operations, accounting for the 



10' 



I 10"' 
c 

3 



10 



XXX' XX 




10 



10 



10 

(OX 



10 



10 



FIG. 1: (Color online) Filter functions for primitive X and 
Xdcg, based on Eq. Q. Upper inset, schematic construction 
of Xdcg [H]- X~ is a X rotation with tt phase shift. Lower 
inset, Amplitude, Fi^x (dotted), and Phase Fi^z (dashed) filter 
functions for the same gates, denoted by color. Tick marks 
same as main panel. 
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FIG. 2: (Color online) Calculated error rates for primitive 
X and Xdcg gates in the presence of noise, (a) Calculate 
error rates in the presence of noise similar to that observed 
in Ref. [5^, appropriate for decoherence due to nuclear spin 
diffusion in solid-state singlet-triplet qubits. Data are plotted 
in dimensionless units of the minimum value, (b) Similar 
to (a), but using a white noise power spectral density up to 
a sharp high-frequency cutoff, (c) Comparison of calculated 
error rates based on analytical filter functions derived herein 
and brute-force numerics. Noise strength in numerics is set 

using A = 0.008 ^ri'^'"'^ , guaranteeing convergence of the 
Magnus expansion for this range of t^. Each data point av- 
erages over 50 randomly generated noise trajectories whose 
statistical properties reproduce S{lu). Error bars are derived 
from the root-mean-square deviation of the individual trajec- 
tory results, but are small compared to the marker size. 



fact that pure dephasing noise can result in both de- 
phasing and polarization damping errors during control 
operations. Included in this class of control sequences 
are dynamically corrected gates, and we have validated 
previous perturbative calculations for DCG performance 
through an intuitive noise-filtering approach. 

We believe that this experimentally accessible and rig- 
orously tested method for understanding the influence of 
classical noise in single-qubit logic operations will prove 
valuable to the community and will open new pathways 
for the development of error-robust quantum control 
strategies. We emphasize that our approach is technol- 
ogy independent and these methods apply to any quan- 
tum system. 
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